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1. Introduction 

In addition to reducing systematic errors in lattice QCD due to finite lattice 
volume, finite lattice spacing and extrapolation from heavy quark masses to 
physical m u and mj, it is important to improve the signal in observables. The 
particular observables we are interested in are matrix elements within hadronic 
states; the corresponding correlation functions from which these are extracted 
are made up of products of quark propagators. Even though at present the 
quark propagators are calculated without specific information on the motion of 
the other valence quarks in a given bound state, it has been shown that the 
signal is significantly improved by using smeared quark sources [1] [2]. 

A number of different smearing methods have been investigated and some 
of these have been reviewed in Ref. [3]. Our choices here are motivated by the 
requirements of a parallel calculation of weak matrix elements using Wilson 
fermions [4]. Our calculation requires two kinds of hadron sources: one that 
produces hadrons with zero momenta and the other that couples to all mo- 
menta. We construct zero momentum hadron correlators using "wall" source 
propagators while the "Wuppertal" source propagators [2] yield hadron corre- 
lators that have overlap with all momenta. In this paper we demonstrate the 
efficacy of these two kinds of correlators by investigating the signal in both the 
amplitudes and the masses obtained from 2-point correlation functions. The 
improvement is exemplified by the results obtained for the spectrum and meson 
decay constants. 

We do, however, find a disconcerting difference between the masses of 
baryons extracted using Wuppertal and wall sources. Even though at this stage 
the difference between the central values is 1-3 a, it has a significant effect on 
the nucleon to rho mass ratio. 

An important prerequisite when using Wilson fermions to calculate matrix 
elements is an understanding of the realization of chiral symmetry. It was 
shown in Refs. [5] and [6] that the conventional continuum current algebra can 
be reproduced on the lattice provided suitably modified chiral Ward identities 
are used. Among other things, this requires multiplying lattice fermion bilinears 
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by renormalization constants to relate them to their continuum counterparts. 
We use the following definitions 



S cont (x) = Z s S(x); S(x) = 4>{x)^{x), 

P cont (x) = ZpP(x); P(x) = lp(x) 75 rP(x), 

V^(x) = Z v V,(x); V,(x) = Mx^^ix), 
A c ° Dt (x) = ZaA^x); A^x) = i>{x) ltll ^{x). 



(1.1) 



If we were working with a chirally invariant regulator and a mass independent 
continuum renormalization scheme, we would have Zy = Za = 1 and Zs = 
Zp = l/Z mass , where Z mass is the mass renormalization. The latter three 
constants have anomalous dimensions and therefore depend on the momentum 
scale. On the lattice with Wilson fermions Zy and Za are finite renormalization 
constants and Zp ^ Zs- Although all these constants can be calculated in 
perturbation theory (see Appendix A), it appears that for present couplings 
(g 2 k, 1) there are large non-perturbative contributions. For Za and Zy , non- 
perturbative calculations are possible, and have been performed by a number of 
authors. These estimates include, in general, large 0(a) effects, and therefore 
different non-perturbative techniques can give different estimates. We use the 
estimates most appropriate for our study, Za = 0.86 and Zy = 0.57 [7]. 

The constants Zs and Zp are more problematic because, unlike Za and Zy , 
they depend upon the choice of continuum renormalization scheme. We need to 
determine these constants if we are to relate the cutoff dependent lattice results 
for the chiral quantities m s and (ipip) m = o to the scale dependent continuum 
values. A convenient continuum scheme is the MS scheme at /j, = 1 GeV. Since 
this is a perturbative scheme, the most natural thing to do is to use the lattice 
perturbative results for Zs and Zp. These are discussed in Appendix A. If per- 
turbation theory is accurate (which is unlikely in the present simulation) then 
one can make unambiguous predictions for m s and (ipip) m = o m the continuum. 

Alternatively, one can proceed in a phenomenological fashion and demand, 
for example, that the lattice results reproduce the value of m s in the MS scheme. 
With this assumption Zp is fixed, and the Ward identities then determine Zs 
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non-perturbatively. The question is then one of consistency: with Zs and Zp 
fixed in this way, can all lattice quantities whose definitions depend upon Zs 
and Zp be made to agree with the MS scheme? In quenched QCD we can 
expect this to be only approximately true. We address these questions further 
in Sees. 7 and 9. For the present we simply state that our data do favor values 
of Zp and Zs very different from their perturbative values. 

This paper is organized as follows: in Sec. 2 we describe the lattices used 
and our method for extracting physical parameters from 2-point correlators. 
The types of quark sources used to construct the hadron correlators are defined 
in Sec. 3. We present an analysis of the spectrum in Sec. 4 and of the meson 
decay constants, f n and fy 1 , in Sees. 5 and 6. The calculation of the current 
quark mass and chiral condensate is presented in Sec. 7, and of the second 
moment of the pion distribution amplitude in Sec. 8. The extrapolation to the 
chiral limit is discussed in Sec. 9, and we end with conclusions in Sec. 10. The 
perturbative results for the renormalization constants are given in Appendix A. 

2. Lattice parameters and fitting procedure 

Our statistical sample consists of 35 lattices of size 16 3 x 40 at (3 = 6.0, 
generated as two independent streams. The first stream consists of 14 lattices, 
generated using a pseudo heatbath algorithm, which are separated by 1000 
sweeps. The second stream consists of 21 lattices generated using a combination 
of over-relaxed and Metropolis algorithms, and are separated by 300 combina- 
tion sweeps. These sets of lattices have been used previously for spectrum and 
weak matrix element analysis using both Wilson and staggered fermions [8]. 

The major thrust of this study is to explore improved numerical techniques 
for Wilson fermions. For this purpose we use only two values of the quark mass, 
k = 0.154 and 0.155. Translated to physical units these correspond to pion mass 
values of roughly 660 MeV and 540 MeV, respectively. The criterion used to 
judge whether the quark propagator has converged is R = | M ^~^ | 2 , where 
X is the solution, <f> is the source vector and M is the Wilson operator. We 
find no significant difference in the long time tail of propagators when reducing 
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R from 10" 15 to 10" 18 , at either quark mass. To be conservative, we adopt 
R = 3 x 10~ 16 as our convergence criterion. 

To extract the amplitudes and masses from the long time behavior of the 
2-point correlators, we make fits assuming that the lowest mass state dominates 
the correlation function. To ensure this we first examine the effective mass plot 
for a plateau and then make a single mass fit selecting the range of the fit 
based on the following criteria: (a) t m ; n always lies in the plateau, (b) t max is 
selected to be as large as possible consistent with a signal. In most cases we 
find that the central value obtained from the fits is the same with and without 
using the full covariance matrix. In some cases we cannot use the full range of 
the plateau because the covariance matrix is close to being singular. In such 
cases the problem is not that we cannot invert the covariance matrix but that 
the result is very sensitive to the range of the fit (the central value can change 
by one or more standard deviations on the addition of a single point to the 
fit range). Our tests using subsets of the data show that this instability is a 
result of inadequate statistics. We, therefore, reduce the range of the fit to get 
a stable result. For example, even though the signal in the pion correlators lasts 
until t pa 40, we use t max < 25 in order to avoid fits based on almost singular 
covariance matrix. 

All error estimates are obtained using a single-elimination jack-knife pro- 
cedure. Previous analysis [8] [9] leads us to believe that the lattices used are 
sufficiently decorrelated for this method to be adequate. 

3. Quark propagators and hadron correlators 

The calculation of quark propagators is done on lattices doubled in the 
time direction, i.e. 16 3 x 40 — ► 16 3 x 80. We use periodic boundary conditions 
in all four directions. These propagators on doubled lattices are identical to a 
linear combination of propagators calculated with periodic (P) and antiperiodic 
(A) boundary conditions on the original 16 3 x 40 lattice. For the source on 
time slice 1, the forward moving solution (time slices 2-40) is F = (P + A)/2 
while the backward moving solution (time slices 80-42) is B = (P — A)/2. We 
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find no numerical difference between calculating F and B directly or from the 
combination of P and A. We use doubled lattices because our weak matrix 
element calculations require both forward and backward moving solutions. We 
will show that doubling also helps in the spectrum analysis. 

We calculate quark propagators using two different types of smeared 
sources: "Wuppertal" and "wall". The wall source for a given spin and color 
consists of a delta function at each point on time slice t = 1, with the wall 
time slice fixed to Coulomb gauge. The hadron correlators built from these 
propagators have overlap with only zero-momentum states. 

The Wuppertal source is the covariant solution of the 3-dimensional Klein- 
Gordon operator K with a delta function source ati = j/ = z = lon time slice 
t = 1 

3 

K(x; y) = 6 XtV - k kg ^ (U (x, i)6 x +i, y + ) , (3-1) 

8 = 1 

where the parameter kkg controls the size of the smearing. The implementation 
of this smearing method is the same as that used in Ref. [10], and the details 
have been presented there. We choose kkg = 0.181, for which the size of 
the smeared source (defined in Ref. [10]) varies between 4.1-4.5. A statistical 
problem when using the Wuppertal source is that, as one increases kkg, the 
configuration to configuration fluctuations in the hadron correlators increase. 
For our value of kkg, the fluctuations at long time separation are less than a 
factor of 5. 

We also construct a third kind of quark propagator which is smeared at 
both the source and sink points. This is obtained by applying the inverse of the 
Klein-Gordon operator to the result of Wuppertal source propagator at each 
time slice. 

From the three kinds of quark propagators we construct three kinds of 
hadron correlators: (a) wall source and local (point) sink (labeled henceforth 
as LW), (b) Wuppertal smeared source and local sink (LS) and (c) Wuppertal 
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smeared source and sink (SS). We use the following notation to define the 
2-point hadron correlators: 

T LW (t) = <O|0 local (*)J wall (O)|O) * ~°° ( Q l 0lOCal |^l Jwa "l°) e -M« ^ 
T LS (t) = <0|O local (*)J smeared (0)|0} * ~°° ( Q l 0lOCal | fe )W smeared l°) e -Mt t 
T ss (t) = (0|O smeared (t)J smeared (0)|0) * ~°° (O|0 smeared |fe)(fe|J smeared |O) e _ M ^ 

W (3-2) 

where \h) is the appropriate hadronic state that saturates the 2-point correlator 
at large t and M its mass. The operator J creates the hadron and is always 
constructed from smeared (Wuppertal or wall) quark sources in our calculation, 
and O is the operator used to destroy it. The projection onto a definite mo- 
mentum state is always done at the sink time slice. We compare the efficacy of 
these correlators for the three types of quark sources in the next section. 
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4. Spectrum 



The flavor non-singlet meson and the baryon interpolating operators we 
have used are as follows: 



7T2 : 
P : 

P2 ■ 
a : 
ai : 

h ■■ 
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7i 



7i 



*72 b\ c 
U Mi 



(4.1) 



where C = 7472 is the charge conjugation matrix. We specify the total angular 
momentum and the parity for the baryons (the opposite parity baryons are 
obtained from the backward moving solutions using the same operators). For 
the A we use two different Lorentz structures (A^ 2+ has Jz = ±3/2 while 
A^ 2+ has Jz = ±1/2) to check for systematic errors. The correlators for 
the two Jz terms are averaged to improve the statistics. The subscript on u c 
specifies the Dirac index of the spinor. The results for the masses in lattice units 
along with the range and xV^df of the fit are given in Tables 1-4. The effective 
mass plots for each of the three sources at k = 0.155 are shown in Figs, la, 
b and c. The corresponding plots at the heavier quark mass, k = 0.154, have 
smaller errors. 
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We find that the effective mass, m e ff(t), converges from above for SS (as it 
must with a positive definite transfer matrix) and for LS correlators, and from 
below for the LW correlator. An agreement between the different methods 
therefore provides a good test of whether the mass estimate is asymptotic. 

In the pion channel we find that m e ff(t) reaches a plateau at t = 10 for 
LS correlators, at t = 8 for SS correlators, while the LW correlators stabilize 
even earlier at t = 6. The onset of the plateau is independent of the quark 
mass in the two cases considered. In the plateau region the three different 
sets of correlators give consistent results. The signal in all three cases extends 
all the way across the lattice and the errors stay roughly constant with t. A 
theoretical analysis of why this occurs for the pion and not for other states has 
been given by Lepage [11]. The errors in m e ff(t) are comparable for LS and 
LW correlators and roughly a factor of 2 larger in SS. In Ref. [10] we gave 
an argument for why, with Wuppertal smearing, the errors in Tss are larger 
than those in T^s- On the other hand, an analysis of the eigenvalues of the 
covariance matrix shows that the time slice to time slice correlations in Tss are 
significantly smaller. The interplay of these two effects makes the errors in the 
fit parameters comparable in the two cases. 

Within the plateau region we have investigated the dependence of the fit 
parameters on the fit range used and on whether or not we use the full covariance 
matrix. We find that while the variation in M w is no greater than 1%, the result 
for the amplitude can change by as much as 12%. This is one reason why the 
extraction of matrix elements is prone to large errors. 

We find that the operator 7T2 has a better overlap with the pion and m e ff(t) 
reaches a plateau earlier. On the other hand the errors in m e ff(t) are much 
larger. A comparison of signals for the two channels can be made from Figs, lb 
and 2a. 

In Fig. 2a we plot m e ff(t) for the four lowest momentum states coupled to 
the A4A4 LS correlator. It is desirable to get a good signal in non-zero momen- 
tum channels as these correlators play an important part in the calculation of 
many matrix elements, for example structure functions and form factors. The 
data show a reasonable plateau for p = (0, 0, 2tt/16) and the possible beginning 
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of one for p = (0, 27r/16, 27r/16). We compare the resulting energy estimates 
with the continuum dispersion relation in Fig. 2b. 

The signal in the rho correlators is also very satisfactory. There is a plateau 
between t = 10 and t = 30 for all three types of correlators, but unlike the case 
of the pion the errors do increase with t. The final estimates from the three 
sets are consistent. We find that the signal is a little better for the p than for 
the p2 operator. 

There is a reasonable signal in the nucleon (N 1 I 2+ ) channel with a plateau 

3/2 + 

extending between t = 10 and 18. We find that the plateau in the A 2 

3/2 + 

channel has a marginally better signal than in the A/ , and for this reason 
we quote results using A2. The big surprise is the systematic difference in the 
baryon spectrum as extracted using the Wuppertal and wall sources. In Fig. 3 
we show an effective mass plot comparing the LS and LW correlators for both 
nucleon and A2. (We emphasize that either Wuppertal correlator, SS or LS, 
could equally well have been used for the comparison with LW .) Both types 
of correlators appear to reach plateaus, though from opposite directions. The 
problem is that the asymptotic values differ, the LW results lying below the LS. 
The same is true for both values of the quark mass. To explore the statistical 
significance of this difference we performed a jacknife analysis of the difference 
of the masses. We find that the difference for the nucleon is only 1 a, while that 
for the A is 2-3 a. Even though the effect is of marginal significance in each 
channel, the fact that in all three baryon channels the LW masses lie below the 
LS and SS masses suggests that the difference is systematic. 

Further evidence for this conclusion comes from a comparison of our results 
with those from the APE collaboration [12], both of which are collected in Ta- 
ble 5. (Our results use the operators tt, p, N and A2.) The APE collaboration 
uses yet another type of source, the "multicube", which is, in a sense, inter- 
mediate between wall and Wuppertal sources. The meson masses are in good 
agreement, while the APE results for baryons lie in between ours for Wuppertal 
and wall sources. 

A similar effect has been observed in hadron mass calculations using stag- 
gered fermions [13] [14]. Meson masses agree, while baryon mass estimates with 
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multicube source again lie a couple of standard deviations higher than those 
with wall sources. 

If this is a systematic effect, it means that the extracted mass estimate is 
not asymptotic. Any such difference has important implications for the M^/M p 
ratio. Calculations over the last two years have found a smaller value for this 
ratio than previous lattice results. This improvement has come, however, from 
calculations that use wall or similar sources [3]. Clearly this is an issue deserving 
further study. 

For the positive parity mesons and the negative parity baryons the signal 
is marginal. In most cases the plateau is noisy and extends only over ~ 5 time 
slices. For this reason we give the results of the fits but do not present a detailed 
analysis of the data. 

The signal in the 7r and p channels extends beyond half the length of the 
original gauge configuration (t = 20), so we can address the question raised by 
the HEMCGC collaboration about fluctuations in the effective mass induced by 
doubling the lattices [15]. The plateaus in our m e ff(t) plots do show fluctuations, 
but these are not correlated with the lattice periodicity. Furthermore, if we 
analyze different subsets of lattices we do not find fluctuations at the same values 
of t. In addition, we have calculated the spectrum using quark propagators 
on the undoubled 16 3 x 40 lattices with periodic boundary conditions in the 
time direction. These results are included in Table 5. There is no significant 
difference between the mass estimates from doubled or undoubled lattices. In 
sum, our data shows no evidence that the fluctuations are anything other than 
statistical. It remains possible, however, that systematic fluctuations would 
appear if the statistical errors were reduced. 

5. Pseudo-scalar decay constant f w 

The definition of the pseudo-scalar decay constant on the lattice is 
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where Za is the axial current renormalization, and we are using the convention 
that the experimental value is f w = 132 MeV. In order to extract f w we study 
two kinds of ratios of correlators: 

T L s(t) t^oo (0|V° ca V) 



x \ / r . . / 4\ i r\\ a smeared | 

(5.2) 



Tss(t) (0\A 4 S 



7T 



R2it) ~ T SS (t) ~ m e ' 

where the T are defined in Eq. (3.2). In the case of R\ we have to extract 
(0|j44 Smeared |7r) separately from the Tss correlator. For each of the two ratios, 
Ri and R2, the smeared source J used to create the pion can be either 7r or 
7T2. Thus we calculate f w in four different ways, which we label as (using 
ratio i?i with J = tt), j\ (using ratio R\ with J = 7^), /| (using ratio R2 with 
J = 7r), and (using ratio i?2 with J = 7^). The results are given in Table 6. 
Errors are estimated by applying single elimination jackknife to the estimate of 
f w itself. We note in passing that the "naive" estimate obtained by combining 
the errors in each of the fit parameters in quadrature results in errors that agree 
with quoted errors to better than 10%. 

We find that all four methods give consistent estimates of f w . Since we 
have to combine different correlators in order to extract f w we select the fit 
range based on the following criteria: (1) goodness of the fit, (2) presence of 
a plateau with a similar mass estimate from each of the individual correlators. 
Otherwise, as stated before, there can be a large variation in the estimate of 
the amplitude. The quality of the signals for the correlator ratios leading to 
and /| are exemplified by Figs. 4a and 4b. 

Bernard tt al. [16] have calculated f w using point sources at the same two 
values of k and using a sub-set of the lattices analyzed by us (they use every 
other lattice). They get f w = 0.094(7) and 0.086(9) at k = 0.154 and 0.155 
respectively, in good agreement with our values. We also find agreement with 
recent point-source results of the QCDPAX collaboration [17]. The results of 
the APE collaboration [12] are in error by a factor Ve" M ' [18]. Aft er correcting 
for this factor their value at k = 0.155 (obtained using extended cube sources) 
is consistent with our corresponding result f£. 
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We have also calculated f w at non-zero momentum, and find that its value 
is stable to the addition of one unit of momentum. The results are included 
in Table 6. The overall quality of the signal is good even though the plateaus 
in the effective mass are much shorter than those for the p = (0, 0, 0) case. At 
p = (0, 27r/16, 27r/16) the errors are larger and the determination is far less 
reliable. 



6. Vector decay constant jy 

We use the local vector current Vi to define the dimensionless number fy 1 

as 

m\P) = ^- (6.1) 

This matrix element can be extracted in two ways analogous to Eq. (5.2). We 
use Vi to both create and annihilate the vector meson. The results are given in 
Table 7, where we use the non-perturbative estimate Zy = 0.57 obtained using 
2-point correlators in an earlier calculation of fy at the same [3 [7] We point 
out that due to the existence of a conserved vector current on the lattice, this 
estimate of Zy is free of 0(g 2 ) ambiguities. 

The quality of the signal in correlator ratios R\ and R2 is very good as 
shown in Figs. 5a and 5b. The final data is shown in Fig. 6 where for com- 
parison we have also included results from the Wuppertal [19], APE [12] and 
QCDPAX collaborations [17]. The experimental points have been taken from 
Ref. [19]. Our results lie significantly below Wuppertal estimates and are in 
agreement with the results obtained by the QCDPAX collaboration and the 



APE collaboration (the latter after correction by a factor \/e M f [18]). The 
data by the QCDPAX collaboration, obtained using point source propagators 



measured on 160 lattices of size 24 x 54, show that smeared and local 
yield consistent results once systematic errors are under control. 
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7. Quark masses and the chiral condensate 

In order to extract the quark masses and the condensate with Wilson 
fermions one must understand how continuum current algebra relations are 
realized on the lattice. This was explained in Refs. [5] and [6], and we recall 
here the relevant results. 

In the continuum, the PCAC relation is 

(ald.A;^) = ( a \^{^-,m q } l5 m, (7-1) 

where \a) and \f3) are physical states, and A a is a flavor Gell-Mann matrix. 
Individually m q and t/'TsV' are scale dependent, and one must choose a particular 
scheme to precisely define them. The standard choice is the MS scheme at a 
scale ii = 1 GeV [20]. 

The lattice relation corresponding to Eq. (7.1) is [6] 

(a\Z A d,A;\p) = Z P (a&{^-,m q } l5 iP\p) +0(a), (7.2) 

where m q , the quark mass, is related to the bare lattice quark mass by 

m q = Z mass m , m =- . (7.3) 

In this and the following equations we assume degenerate quarks. We stress 
that Z mass is the mass renormalization and differs from the finite constant Z m 
defined in Ref. [6] (Z m = ZpZ mass in our notation). In Ref. [21] it was shown 
that for m q — ► 0, Z mass = Zg 1 . In Appendix A we rederive this result in 
perturbation theory to all orders, and argue that it remains true for general 
m q in the continuum limit. For this reason, we will express all our results in 
terms of Zs and Zp by requiring -Zmass = Zg 1 identically. It should be borne in 
mind, however, that for non-zero lattice spacing there are O(ma) corrections to 
the individual Z factors which depend on the initial and final states. Thus Zs 
factors extracted from different matrix elements may differ by terms of O(ma). 

For Wilson fermions, the absence of chiral symmetry means that Zs and 
Zp differ, and as discussed in the introduction they must either be calculated 
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perturbatively or else fixed by some phenomenological requirement, since the 
continuum quantities we are comparing to are defined perturbatively. Their 
ratio Zp/Zs, however, is finite and can be extracted non-perturbatively, as we 
now describe. 

From Eq. (7.2) we can extract m q Zp j Z a ■ To do this we calculate the 

_M JL (0|A 4 (f)J(0)|0) t^co Zp_ a 
2 <0|P(f)J(0)|0> ~ Z A m * 

1 (0|d4^4ft)J(0)|0) t^oo Zp_ c ( ' ) 

2 (0|P(*)J(0)|0) ~ Z A m «' 

where P and A 4 are the local operators and the smeared source J = ir. Thus, 
given Za and assuming m q = Z^ 1 m 1 ? att , we can make a non-perturbative evalu- 
ation of Zp /Zs- Two other estimates, m b q and m d q , are obtained by substituting 
J = TT2- All four methods give consistent results, as shown in Table 8. The 
quality of the signal is displayed in Figs. 7a and 7b. Also given in Table 8 are 
the values for m' g att and the results for Zp/Zs- Our values are consistent with 
the earlier estimate of Ref. [7], Zp/Zs ~ 0.7. 

To calculate the chiral condensate we use two variants of the method sug- 
gested in Ref. [6]. This is based upon the continuum Ward Identity 



(#) WI = (015(0)10) = lim m q / d 4 x(0\P(x)P(0)\0) , (7.5) 

where (t/>t/>) WI is the chiral condensate per light flavor. The lattice equivalent 
of this is [6] 

(^) WI = lim m q y"(0\ZpP(x)Z P P(0)\0). (7.6) 

m q — ^0 ^ — d 

x 

Using Eq. (7.4) one can rewrite this as 

^-^^^^^^ <"> 

where the ratio of correlators is evaluated at large t. We cannot use (7.7) since 
we only have LS and SS correlators available, so instead we use 

(H^_ 1 {o\d 4 A 4 (t) j (o)\o) 



ZpZ A 



m™0 2 (0|P— ed^ J(0) |0) 2J°I^ > P (°)l°>- ( ? - 8 ) 
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This is equivalent to Eq. (7.7) if the pion pole dominates the sum on the right 
hand side, which occurs in the limit m q — ► 0. The results are given in Table 9 
for the two choices of the source J. A typical example of the quality of the 
signal is shown in Fig. 8a for J = ir. 

A variant of this method is to assume pion dominance of Eq. (7.7), and 
derive a lattice version of the Gell-Mann, Oakes, Renner relation [22] [20] 

(#} GMOR = Km Jl^l. (7.9) 

x ' m q ^o Am q v ' 

Using Eqs. (5.2) and (7.4), this can be extracted from the combination of cor- 
relators 

(0|P(t)J(0)|0)(0|A 4 (t)J(0)|0) t _o W) GMOR Mri 

(o|j(*)j(o)|o) ~ z P z A ' { ■ ' 

where the pion source can again be J = 7r or 7T2. (To derive this result one must 
bear in mind the normalization factors of 2M W , as shown explicitly in Eq. (5.2).) 
The results using this method are also given in Table 9, and in Fig. 8b we show 
a typical m e ff(t) plot for the above ratio. 

These two methods should only agree for m q = 0, in which limit they give 
the condensate. In an expansion 

{#) = (#) m ,=o + «n f + -- -, ( 7 -H) 

the linear (and higher order) terms are contaminated by lattice artifacts. We 
also include in the Table 9 the results of linearly extrapolating (ipip) to k c = 
0.15704. It is very encouraging to see that the four values agree with each other 
within our statistical resolution, and this leads us to believe that the results 
are physical. Therefore we compare them with those obtained using staggered 
fermions and with the experimental values in Sec. 9. 
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8. Second moment of the quark distribution amplitude in the pion 

The quark distribution amplitude of a hadron is a wave function describing 
the distribution of the hadron momentum between valence quarks near the light 
cone limit. For the pion, the second moment of the distribution amplitude, (£ 2 ), 
parameterizes the matrix element of an axial operator with two derivatives: 

(0\A^ p (0)\ir(p)) = Z^f^PuPuPp, 

(8.1) 

A [lvp {x) = (-i) 2 ip(x)j 5 j^D v D p ip(x) - traces. 

Eqn. (8.1) holds independent of the symmetrization of Lorentz indices, though 
for certain analyses such a symmetrization is desirable in order to project out 
an operator of definite twist. On the lattice, the renormalization constant Z 
is unknown, and depends on the Lorentz indices. In one loop perturbation 
theory, assuming that the dominant contribution to the renormalization comes 
from the tadpole terms which are independent of the Lorentz indices, one can 
estimate Z ~ 1.3. There is no reason to trust this perturbative result, however, 
and the results given below support a value substantially larger. 

Using the lattice transcription of operators given in Ref. [23], we have 
measured ^ 2 ) latt = Z _1 (^ 2 ) from the following correlator ratios: 



R [433] _ (0|^[433](P,<)J(0)|0) _ ,, 2 . latt n n 

Kj ~ {0\A 4 (p,t)J(0)\0) " (4 > P3P3 > 



ff [343] _ (0|A343](^)J(0)|0) _ 2 latt . 

Rj ~ (0\A 4 (p,t)J(0)\0) " (0 P3P3 > (8 ' 2) 



r>(433) 

Kj ~ {0\A 4 (p,t)J(0)\0) " (4 > P3P3 ' 



where J = 7r, 7T2 and p = (0, 0, 27r/16). Square brackets around indices in- 
dicate that an appropriate combination of operators has been taken to non- 
perturbatively subtract quadratic divergences [24] and round brackets indicate 
that symmetrization over the Lorentz indices has also been performed: 

^4[433] = ^4-433 — ^4.411, 

^4[343] = ^343 — ^4.141, (8-3) 
^4(433) = (/4[433] + 2A[ 343 ])/3. 
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Our results are given in Table 10. Examples of the quality of the data are shown 
in Figs. 9a and 9b. In general using 7r as the pion source seems to give better 
results: with 7T2 it takes longer to reach the asymptotic plateau. 

For the symmetrized operator, j4( 433 ), the results are consistent with a pre- 
vious analysis [23] in which configurations with dynamical quarks with masses 
in the range m s < m q < 3m s were used. This suggests that (,ff 2 ) latt & 0.1. 
However, examination of the two unsymmetrized operators separately reveals 
a surprising feature: (,ff 2 ) latt from the operator with a time derivative, j4[ 343 ], is 
much larger than from that with only spatial derivatives, j4[ 433 j. 

There are various possible explanations of this discrepancy, but note in 
particular that j4[ 433 ] and j4[ 343 ] renormalize differently (they do, of course, 
mix). This can resolve the problem only if there are large non-perturbative 
contributions to the renormalization constants. Lattice spacing corrections may 
also play a large role, and repeating these calculations with an improved fermion 
action should be very helpful. 

We have also tried to determine (^ 2 ) latt from A423 at p = (0, 27r/16, 27r/16). 
The results are given in Table 11, though we are less confident of the signal at 
this higher momentum. We associate the poor signal with the lack of a plateau 
in the effective mass plot for the pion correlator with two units of momentum. 
The renormalization of j4 4 2 3 is independent of the ordering of the indices, and 
the data do not show any significant difference. 

9. Lattice scale and comparison of results 

In order to study the chiral limit, we extrapolate our results for the various 
quantities to the value of k at which the pion mass vanishes. This is shown in 
Fig. 10, where we have made separate extrapolation for baryon masses obtained 
using Wuppertal (superscript a) and wall (superscript b) sources. We extrapo- 
late using m' 3 att . We find k c = 0.15704, in agreement with the result given in 
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Ref. [12]. The results of these fits are: 



Mia 2 


= 0.000(16) - 


f 2.1(3) m 


fir a 


= 0.070(14)- 


f 0.38(27) m 


M p a 


= 0.313(33) - 


f 2.3(6) m 


M%a 


= 0.521(53) - 


f 3.4(1.0) m 


M h N a 


= 0.486(91)- 


f 3.5(1.6) m 


M^a 


= 0.700(72)- 


f 1.9(1.4) m 


M b A a 


= 0.58(13) + 


2.9(2.1) m 



(9.1) 



from which we can extract the lattice scale. Using the experimental values for 
f w , M p , Mjv and Ma we get: 

a -1 CM = 1-9(4) GeV 
a -\M p ) = 2.5(3) GeV 
a- 1 (Ml) = 1.8(2) GeV 

(9.2) 

a~ 1 {M b N ) = 1.93(36) GeV 
a- 1 {Ml) = 1.76(18) GeV 
a~ 1 {M b A ) = 2.1(5) GeV. 

The large errors in the extrapolated values reflect the fact that we have data at 
just two values of k. The a -1 extracted from M p is significantly larger. Also, 
the wall source estimates for the nucleon mass give a slightly better value for 
the ratio M N /M p . 

A linear extrapolation of (ipip) to the chiral limit gives, 

{H) mq =o = -Z P Z A 0.0058(18) a" 3 

(9.3) 

~ -Z P Z A 0.034 GeV 3 , 

where we have used a -1 = 1.8 GeV. Assuming ZpZ A ~ 1, we note that this 
result is approximately three times larger than the continuum value (uu) = 
-0.0114 GeV 3 [20]. To make a meaningful comparison we need to estimate Zp. 
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The 1-loop perturbative result for Zp using Wilson fermions is discussed 
in Appendix A. It is 

Zp = i + f^T ( 6 lQ gM- 22.596), (9.4) 

where C'f = 4/3. This estimate does not address possible 0(a) and other 
non-perturbative effects. We use the continuum MS renormalization scheme at 
l_t = 1 GeV so that Zp = 1 — .22 1 <? 2 for a -1 = 1.8 GeV. This leaves one source of 
uncertainty in applying this formula, namely the choice of g 2 . Using g 2 2, as 
advocated by Lepage and Mackenzie [25], we get ZaZp 0.5. Given the many 
uncertainties in this estimate, what we can say with some confidence is that 
including these corrections moves the lattice estimate closer to the continuum 
result. 

Using Eq. (9.1) we can also attempt to extract the mass of the strange 
quark. Our estimate varies by a factor of two depending on the method we use. 
For example, by demanding that the ratio 2M^-/f 2 attain its physical value 
at m = m s we get m\ att a 0.066, while the ratio M^»/2M^ gives m\ att a 
0.036. Translated into physical units these correspond to Zsm s = 81 MeV and 
44 MeV respectively. If we assume Zs ~ I, then these values are a factor of 
2-3 smaller than the conventional estimate. This discrepancy is similar to the 
results obtained with staggered fermions [13]. Once again, a value of Zs < 1 
would increase these estimates. 

The combination m s (ipip) m = o can be extracted with less ambiguity be- 
cause it involves the finite ratio Zp/Zs, for which we have a non-perturbative 
estimate. Taking Za = 0.86 and Zp/Zs = 0.68, and using the two extreme 
values for m s quoted above, we find —0.0024 GeV 4 and —0.0013 GeV 4 , to be 
compared to the experimental value of approximately —0.0017 GeV 4 . 

10. Conclusions 

We show that both Wuppertal and wall quark sources yield very good 
signals for the 7r and p mesons at p = (0,0,0). There is an unambiguous plateau 
in the effective mass plots for all three types of correlators studied and the mass 
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estimates are consistent. Since the estimates converge from opposite directions 
with the two kinds of sources, consistency of the results implies that we have 
extracted the asymptotic value. We find that the signal in the Wuppertal source 
pion correlator with p = (0, 0, 27r/16) is good enough to allow the calculation 
of matrix elements with non-zero momentum flow. 

There are strong indications that we have not extracted the asymptotic 
value of the baryon masses. The disturbing difference between the results from 
the Wuppertal and wall sources needs to be understood before we can quote 
the baryon masses with confidence. 

The signal in the ratios of Wuppertal source correlators used to extract f w , 
fy 1 , m q and (ipip) is very good. We show how to extract the current quark 
mass and the chiral condensate using smeared sources, although comparing with 
experiment is difficult in the absence of reliable values for Zp and Zs ■ We show 
internal consistency of the results by using different hadronic operators and by 
using different combinations of correlators. We find that our estimate of fy 1 is 
in good agreement with experimental values. 

Lastly, we have calculated the second moment of the pion distribution am- 
plitude, (£ 2 y att . We find a significant difference between the results extracted 
from j4[433] and j4[ 34 3] at p = (0, 0, 27r/16). The average value of (^ 2 ) latt ~ 0.1 is 
consistent with earlier results obtained using lattices with two dynamical flavors 
of Wilson fermions. 
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Appendix A. Perturbative evaluation of Zs, Zp and Z mass 

The perturbative calculations required to evaluate the constants Zs, Zp 
and -Zmass with Wilson fermions have been performed by a variety of authors [26] 
[27] [28] [29] [30] [31]. It is, however, rather difficult to extract the precise values 
from these papers, and for this reason we collect some of the relevant results 
here. We restrict ourselves to r = 1, and take MS as our continuum scheme. 

The calculation of the renormalization constants requires the evaluation of 
the fermion self-energy (wave function and mass renormalization) and vertex 
corrections. The continuum fermion self-energy in the MS scheme has the form 

S^(p,m) = if -pS;p(p,m) + mSf (p,m), (A.l) 

with the 1-loop results 

rl ' xm 2 + x(l — x)" 2 - 



EHp )ra) = A{l + 2 1 dx { l-x)lo g ( Xm +X j;r X)P )}, 
EFV m) = A{2 + 4 jf 1 dx log ( W + ; ( 2 1 -^ 2 )}, 
where A = (/ 2 CV/(167r 2 ) and Cf = 4/3. The lattice self-energy is 



(A.2) 



E latt (p, m) = i S' att + ij -p^ip, m) + m £ 2 att (p, m), (A.3) 

where at 1-loop 
E 1 *" = A(-51.435), 

S'i att (p, m) = A|l3.852 + 2 ^ rfa; (1 - x) log (a 2 (;cm 2 + a;(l - x)p 2 )) } 

+ 0(ma,pa), (A. 4) 

£ 2 att (p, m) = A|l.901 + 4 J dx log (a 2 (a;m 2 + x(l - x)p 2 )) } 
+ 0(ma, pa). 

We have taken the most precise values available for the finite constants [29]. 
These are obtained by numerical integration and are accurate to better than 
one part m the last decimal place. The linearly divergent piece, X;Q att , shifts 
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the position of k c , and plays no role in the following discussion. Neglecting this 
term, and adopting a slightly modified version of the notation of Ref. [27], we 
write the difference between continuum and lattice schemes (in the limit a — ► 0) 
as 

S MS (p,m) - S latt (p,m) = ij -pA Sl + mA S2 . (A. 5) 

At 1-loop, let us define A; = AA^ 1 - 1 . Then 

A^ = -2 \og(ap) - 12.852, A^ = -8 \og(ap) + 0.099. (A. 6) 

With these results the 1-loop mass renormalization constant is 

^mass = 1 + ^(Aj^ - A^), ^ 

= 1 + A(- 6 log( a/ u) + 12.951). 

In Ref. [26] the result for Z mass is in error, because the expression for A^ is 
incorrect. This error was pointed out in Ref. [28], where a result consistent with 
Eq. (A. 7) is obtained. 

Notice that like the A's, the perturbative Z's are defined neglecting terms 
of 0(a) and, in particular, terms of O(ma). 

To calculate Zs and Zp one also needs the difference between the MS and 
lattice vertex corrections for insertions oiipip and i/'TsV'- Following Ref. [27], we 
write these differences as Ai and A 76 respectively. The above mentioned error 
in Ref. [26] propagates into an error in the 1-loop result for Ai = AAj 1 - 1 given 
in Ref. [27]. The result for A^- 1 is correct. More precise 1-loop values may be 
deduced by combining the results of Refs. [31], [30] and [27]: 

A^ = +8 \og{an) - 0.100, = +8 log(a^) - 9.744. (A. 8) 

To extract these results one needs to know that in the dimensional reduction 
scheme used in Ref. [30], both A^ 1 - 1 and A^- 1 are larger by +2 than in the MS 
scheme. 
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With these results at hand we can now calculate the 1-loop renormalization 
constants 

Z s = l + A(A^ + Ag), 

= l + A( + 61og( a/ u)- 
Z P = l + \(AW + A ( i}), 

= l + A( + 61og(a//)- 

Eqs. (A. 7) and (A. 9) imply that, because A^ 1 " 1 = — A^j within the numerical 
errors, the relation Z mass Zs = 1 is satisfied. As mentioned in the text, this 
relation is expected on the basis of the functional integral derivation given in 
Ref. [21]. In fact, as we now discuss, it is straightforward to understand this 
result in perturbation theory. 

To establish the equality Ai = — As 2 at 1-loop we begin by noting that to 
this order 

AA£ = ^ (^(p, m) - ^ tt (p, mj) \ m=0 , (A.10) 

which follows from Eqs. (A. 5) and (A. 6). The crucial observation is that the 
derivative with respect to m inserts minus the scalar density on the internal 
fermion line, and so gives the difference of continuum and lattice vertex dia- 
grams. These are precisely the diagrams which, in the limit m —> 0, define 
— A^ 1 - 1 . Note that terms suppressed by powers of ma (which are present in 
S latt ) are irrelevant since the A's are defined dropping all terms of 0(a). 

This derivation extends straighforwardly to all orders. Taking the deriva- 
tive with respect to m in a general self energy diagram is equivalent to the 
insertion of minus the scalar density in all possible ways, and so by definition 
gives the corresponding set of graphs for — Ai. Note that this is true for any 
r. The crucial assumption is that the difference of continuum and lattice self- 
energies which defines the A's is infra-red finite, and thus has a smooth chiral 
limit. This forbids contributions to the difference of S's which would make the 
derivative infrared singular, i.e. terms of the form log(ma), log(m 2 /p 2 ), p 2 /m 2 , 
etc. Dimensional analysis then shows that the only dependence on m comes 
from polynomials in ma and malog(ma), which give no contribution in the 
continuum limit. 



12.952); 
22.596). 



(A.9) 
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This discussion makes clear that the relation ZsZ mass = 1, which we use in 
the text, holds exactly in perturbation theory. The individual renormalization 
factors we observe in a simulation at finite lattice spacing, however, also contain 
O(ma) corrections, the precise values of which depend upon the external states. 
Consequently, the replacement Z m3SS — ► Zg 1 may have corrections of O(ma) if 
the two are extracted from different matrix elements. 

Finally, we give the finite constant 

Zm ass z P = i + a(aW + a£), 

(A. 11) 

= 1 + A(-9.644), 

which is denoted Z m in Ref. [6]. Note that in order to reproduce our non- 
perturbative result using this formula, one would have to take g 2 k, 4. 
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0.759(31) 


0.708(37) 
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0.361(4) 


0.365(6) 


0.463(6) 


0.460(9) 


0.727(30) 


0.676(34) 


0.714(31) 






0.708(18) 


0.935(23) 


0.767(21) 


0.761(20) 


0.958(25) 


1.009(42) 



Table 1. Meson masses at k = 0.154. For each operator and correlator we give Table 2. Baryon masses at k = 0.154. For each operator and correlator we 

the x 2 P er degree of freedom, the fit range and the mass estimate. give the \ 2 P er degree of freedom, the fit range and the mass estimate. 
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0.910(21) 
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Table 3. Meson masses at k = 0.155. For each operator and correlator we give Table 4. Baryon masses at k = 0.155. For each operator and correlator we 

the x 2 P er degree of freedom, the fit range and the mass estimate. give the \ 2 P er degree of freedom, the fit range and the mass estimate. 



K 


Lattice 


Nconf 


7T 


P 


N 


A 




0.154 


16 3 x 40 


35 


0.365(4) 


0.465(7) 


0.736(11) 


0.82(2) 


0.093(3) 


0.154 


16 3 x 80 


35 a 


0.364(6) 


0.460(7) 


0.737(14) 


0.82(2) 


0.091(4) 


0.154 


16 3 x 80 


35* 


0.361(4) 


0.463(6) 


0.708(18) 


0.76(2) 




0.154 


18 3 x 32 


104 


0.361(1) 


0.463(3) 


0.721(7) 


0.782(10) 




0.155 


16 3 x 40 


35 


0.301(6) 


0.420(12) 


0.663(15) 


0.78(2) 


0.086(3) 


0.155 


16 3 x 80 


35 a 


0.297(9) 


0.411(10) 


0.665(15) 


0.78(2) 


0.087(4) 


0.155 


16 3 x 80 


35* 


0.295(5) 


0.420(10) 


0.634(28) 


0.70(4) 




0.155 


18 3 x 32 


104 


0.295(1) 


0.422(4) 


0.651(10) 


0.723(14) 




0.155 


24 3 x 32 


78 


0.297(3) 


0.428(4) 


0.647(6) 


0.745(15) 


0.083(3) 



Table 5. Quenched Wilson Fermion Spectrum at [3 = 6.0. We present results 
from Wuppertal (denoted by superscript a in column 3) and wall sources (su- 
perscript 6) on 16 3 x 80 lattices separately. Our Wuppertal source results on 
both doubled and un-doubled lattices are the mean of the LS and SS values. 
We include results on 18 3 x 32 and on 24 3 x 32 lattices from the APE collabo- 
ration [12] for comparison. To get the value for f w in lattice units we have used 
Z A = 0.86. 





f a 


J 7T 


J 7T 


J 7T 


k = 0.1540 
p = (0,0,0) 


0.95 
15-22 
0.090(4) 


1.56 
12 - 22 
0.090(7) 


1.33 
12 - 22 
0.092(5) 


1.00 
12 - 22 
0.093(7) 


k = 0.1540 
P = (0,0,1) 


0.73 
10- 16 
0.093(12) 


1.84 
8-16 
0.094(5) 


0.21 

9-16 
0.094(7) 


0.32 
10 - 16 
0.094(8) 


k = 0.1550 
P = (0,0,0) 


1.93 

8-24 
0.086(4) 


1.18 

10 - 20 
0.079(8) 


1.49 

8-28 
0.088(5) 


1.85 
9-21 
0.087(7) 


k = 0.1550 

P = (o,o,i) 


0.92 
8-14 
0.090(7) 


1.23 

8-14 
0.085(4) 


0.93 
7-14 
0.084(9) 


0.23 

9-16 
0.088(8) 



Table 6. Results for the pseudoscalar decay constant, f w , calculated in the 
four ways described in the text, and using Za = 0.86. For each measurement 
we give the \ 2 P er degree of freedom, the fit range and the estimate. 



K 






0.1540 


1.03 

8-18 
0.24(2) 


1.53 
10 - 20 
0.25(2) 


0.1550 


1.19 

8-18 
0.26(2) 


2.84 
10 - 22 
0.26(2) 



Table 7. The value of the vector meson decay constant, fy , calculated in the 
two ways described in the text, and using Zy = 0.57. For each measurement 
we give the \ 2 P er degree of freedom, the fit range and the estimate. 



K 


Z P m a q 


Z P m\ 


Z P m\ 


Z P m d q 


m latt 


Zp/Zs 


0.1540 


1.29 
12 - 25 
0.042(1) 


0.79 
12 - 22 
0.041(2) 


1.76 
10 - 30 
0.043(1) 


0.79 
12 - 28 
0.042(1) 


0.062 


0.68(1) 


0.1550 


1.41 

12 - 20 
0.027(1) 


0.85 
12 - 22 
0.027(2) 


1.52 
12 - 28 
0.028(1) 


1.06 
12 - 28 
0.028(1) 


0.042 


0.66(2) 



Table 8. The value of Zpm q , calculated in the four ways described in the text, 
using Za = 0.86. For each measurement we give the \ 2 P er degree of freedom, 
the fit range and the estimate. Also given are the lattice quark mass and the 
renormalization constant Zp/Zs- 



K 


{4>4>)T 


(^)Z 1 


(#)™ OR 




0.1540 


1.12 

12 - 22 
-0.0080(8) 


0.84 
8-18 
-0.0079(8) 


1.78 
10 - 20 
-0.0147(8) 


1.22 

8-18 
-0.0149(8) 


0.1550 


1.89 
8-22 
-0.0073(7) 


1.48 
10 - 22 
-0.0072(13) 


1.82 
8-18 
-0.0118(7) 


1.07 

8-18 
-0.0118(8) 


0.15704 


-0.0059(26) 


-0.0058(42) 


-0.0060(26) 


-0.0056(29) 



Table 9. The value of the chiral condensate on the lattice. For each measure- 
ment we give the \ 2 P er degree of freedom, the fit range and the estimate. We 
have to multiply all numbers by ZpZa in order to get physical values. The 
result at k c = 0.15704 is obtained by linear extrapolation. 



K 


{a [433] 


<a? 43] 


{a (433) 


<ar ] 


<ar ] 




0.1540 


0.61 
5-10 
0.06(2) 


0.51 
5-10 
0.13(2) 


0.82 
5-10 
0.10(2) 


0.47 
10 - 16 
0.07(4) 


1.27 
10 - 20 
0.16(3) 


1.32 
6-16 
0.11(1) 


0.1550 


0.66 
5-10 
0.06(3) 


0.54 
6-12 
0.12(3) 


0.91 
5-12 
0.10(2) 


0.16 
10 - 15 
0.09(5) 


1.70 
7-15 
0.18(2) 


1.37 
7-15 
0.11(2) 



Table 10. The value of (^ 2 ) latt from A 433 calculated in the six ways described 
in the text. For each measurement we give the \ 2 P er degree of freedom, the 
fit range used and the estimate. 



K 


<a 4 . 23 


(e)r 


<a^ 423) 


<a 4 r 


(erif 


<ar } 


0.1540 


1.34 
6-12 
0.14(2) 


0.70 
6-12 
0.10(2) 


1.62 
6-15 
0.12(2) 


2.09 
8-12 
0.08(4) 


2.19 
8-15 
0.11(3) 


3.58 
7-14 
0.11(2) 


0.1550 


0.82 
6-12 
0.15(4) 


1.17 
5-10 
0.13(3) 


1.05 
6-12 
0.12(3) 


1.77 
7-12 
0.08(4) 


2.02 
8-15 
0.09(3) 


3.66 
6-12 
0.15(3) 



Table 11. The value of (,ff 2 ) latt from A423 calculated in the six ways described 
in the text. For each measurement we give the \ 2 P er degree of freedom, the 
fit range used and the estimate. 
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Fig. la. The effective mass plots for the ir, p, nucleon and the A at k = 0. 
using SS correlators. 
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lb. Same as in Fig. la except the data are for LS correlators, 
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lc. Same as in Fig. la except the data are for LW correlators. 
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Fig. 2a. The m e ff(t) plot for 1:2^2 LS correlator for the 4 values of momentum, 
p= (0, 0, 0), (0, 0, 1), (0, 1, 1), and (1, 1, 1), at k = 0.155. 
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Fig. 2b. Comparison of lattice pion spectrum with the continuum dispersion 
relation E 2 = m 2 + p 2 . Calculations were done only for the four lowest values 
of momentum. We do not quote a result at k = 0.155 for p = (1, 1, 1) as there 
is no credible plateau in m e ff(t) (see Fig. 2a). 
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Fig. 3. Comparison of the m e ff(t) plots for the baryons at k = 0.154 using 
LW and LS correlators. 
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Fig. 4a. The ratio of correlators R\ at k = 0.154 from which /°, as defined 
in the text, is extracted. 




Fig. 4b. The m e ff(t) plot for the ratio of correlators R2 at k = 0.154 from 
which /£, as defined in the text, is extracted. 
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Fig. 5a. The ratio of correlators R\ at k = 0.154 from which 1/ fy, as defined 
in the text, is extracted. 
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Fig. 5b. The m e ff(t) plot for the ratio of correlators ^ at k = 0.154 from 
which l/fy, as defined in the text, is extracted. 
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Fig. 6. Comparison of the lattice estimates of l/fv with experimental data. 
Points labeled by + are the corrected data from Ref. [12], o from Ref. [17], x 
from Ref. [19] and o from this calculation. Experimental points are indicated 
by □. 
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7a. Plot of the ratio of correlators at k = 0.155 from which m b q , 
in the text, is extracted. 
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7b. Plot of the ratio of correlators at k = 0.155 from which m d q , 
in the text, is extracted. 
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Fig. 8a. The plot for the ratio of correlators yielding (t/>t/>) WI at k = 0.154. 
The result has to be multiplied by 4k 2 to take into account our normalization 
of quark propagators. 



0.5 



0.4 



0.3 



0.2 



~i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — r 
<P A 4 > LS <A 4 A 4 > LS /<A 4 A 4 > ss 



T 



x 



X ' "X"-Mr 
X *X-^ 



X )( 



•x 



X 



XX- 



X 



X, 



X 



" X 
Xx- 



>c 



X 



X 



-X 



X 



\ I I I I I I I I I I I I I I I I I I I L 

10 20 30 

t 



40 



. 8b. The m e ff(t) plot for the ratio of correlators yielding (t/>t/>) GMOR at 
0.155. 
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Fig. 9a. Plot of the ratio of correlators, R w , at k = 0.155 from which the 
lattice value of (^ 2 ) is extracted. 
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Fig. 9b. Plot of the ratio of correlators, R W2 , at k = 0.155 from which the 
lattice value of (^ 2 ) is extracted. 




Fig. 10. Linear extrapolation of the lattice masses to the chiral limit. We 
show both sets of baryon masses, those obtained using Wuppertal (superscript 
a) and using wall (superscript 6) sources. 



